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Introduction 

A new hypersonic test facility which can simulate conditions 
typical of atmospheric flight at Mach numbers up to 2 0 is currently 
under study at the NASA/LaRC Hypersonic Propulsion Branch. In the 
proposed research it was suggested that a combustion augmented 
electrothermal wind tunnel concept may be applied to the planned 
hypersonic testing facility. The purpose of the current 
investigation is to evaluate some candidate working fluid 
formulations which may be used in the chemical-electrothermal wind 
tunnel . 

Work Done 

The efforts in the initial phase of this research were 
concentrated on acquiring the code used by GASL to model the 
electrothermal wind tunnel and testing it using the conditions of 
GASL simulation. The early version of the general chemical 
kinetics code (GCKP84) was obtained from NASA and the latest 
updated version of the code (LSENS) was obtained from the author 
Dr. Bittker. Both codes are installed on a personal computer with 
a 486 25 MHz processor and 16 Mbyte RAM. Since the available 
memory was not sufficient to debug LSENS, for the current work 
GCKP84 was used. 

The Effect of N0 2 Reactions: As a first step we tried to reproduce 
the results obtained by GASL for air-argon mixtures. For this 
purpose the GASL kinetic and thermodynamic data provided in 
reference (1) are used. Although kinetic constants for reactions 
involving N0 2 (reactions 6 and 7 in the GASL table) were given in 
Table 7 of reference (1), it was indicated in the report that these 
reactions were not considered in the computations. Private 
communication with O. F. Rizkalla confirmed this and the reason 
given for the exclusion of these reactions were the unavailability 
of high temperature thermodynamic data and the expectation that 
their influence on the test gas composition will be negligible. 

We have obtained high temperature thermodynamic data for N0 2 
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and several other species of importance to our analysis 2 . The data 
obtained contains two sets of coefficients for the ranges 200 to 
1000 K and 1000 to 10,000 K. For some species the upper limit is 
20,000 K. Since GASL has used three temperature interval fit for 
their thermodynamic data, the N0 2 data could not be incorporated to 
their thermodynamic data, and the McBride data was used in our 
studies as a complete set. 

The initial conditions for the reactive nozzle flow were taken 
from Table 6 of reference (1). A 5° nozzle with the specified 
throat area was used. The nozzle area, A in cm 2 , is given by 

A = 0.0240466 X 2 + 0.198988 X + 0.41166 


In our simulations the reactions and accompanying rate data 
given in Table 7 of the GASL study 1 are used. For some runs the 
kinetic data for reaction 6 was replaced by the data from a test 
case supplied by Dr. Bittker with LSENS 3 , and the results shown in 
the second row were obtained using the high temperature kinetic 
data for reactions 6 and 7 reported by Jachimowski^ . The result of 
our simulations are compared with the GASL results 1 in Table 1. 

Considering the fact that the GASL results are read from a 
graph, the NO, O, and 0 2 mole fractions have been reproduced by 
using both the GASL thermodynamic data and the McBride data. In 
fact, McBride Thermodynamic data appears to produce results closer 
to those of GASL. The only major discrepancy between our 
simulations and the GASL results is the nozzle length required to 
produce an exit Mach number of 16. The nozzle length from Figure 
16 of reference 1 is about 8.9 ft while our simulations indicate a 
nozzle length of about 8.7 m. 

The first 6 rows of the table show the effect of N0 2 kinetics 
on the composition of the gas in the nozzle. The main effect of 
the N0 2 reactions is on the oxygen atom concentration. When these 
reactions are excluded from the calculations the oxygen atom mole 
fraction at the nozzle exit is 1.24E-2 while with the N0 2 reactions 
the oxygen atom mole fraction falls below 7.25E-5 for all rate 
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data. The effect of N0 2 reactions on the NO mole fraction is not 
that significant. The NO mole fraction is 4.84E-2 when the 
reactions are excluded and about 5 . 6E-2 when they are included. 


Table 1. Effect of NO, reactions on the test gas composition 

with air as the working fluid. 


Reaction 

Information 

Conditions at 800 cm (y i in mole frac.) 

T 

(K) 

y N0 

y 0 

y 02 

M a 

N0 2 Reactions Included 
GASL Reaction Data 

271.66 

5.592E-2 

5.625E-7 

0.1813 


N0 2 Reactions Included 
Data from Jachimowski 

270.90 

5.604E-2 

7.253E-5 

0.1813 


N0 2 Reactions Included 
Bittker's Data for 6 

267.04 

5.648E-2 

6 . 883E-14 

0.1811 

15.477 

Reaction 7 Excluded 
GASL Data for 6 

322.81 

3.653E-2 

2.715E-3 

0.1782 


Reaction 7 Excluded 
Bittker's Data for 6 j 

291.57 


1.387E-3 

0.1777 


N0 2 Reactions Excluded 
McBride Thermo Data 

263.20 

4.838E-2 

1.246E-2 

0.1778 


N0 2 Reactions Excluded 
GASL Thermo Data 

260.8 

5 . 106E-2 

1.067E-2 

0.1776 

WM 

GASL Results at Nozzle 
Exit 


bo 

1.9 E-2 



N0 2 Reactions excluded 
McBride Data-Exit Cond. 

247.44 

4.838E-2 

1.244E-2 




Simulation Conditions and Procedure: Results have been obtained at 
Mach 16 free jet conditions for air, air containing 10% ammonia, 
air containing water equivalent to the hydrogen in 10% ammonia 
mixture, and a mixture obtained by replacing some nitrogen in air 
by ammonia to give 10% ammonia. Computations for air and air + 10% 
NH 3 mixtures were done for 5 and 10 degree nozzles while for the 
others only a 5 degree nozzle was used. The total enthalpy was 
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3159.5 cal/g and the mass flow rate was 13302 g/s for all 
computations. The specified total enthalpy is assumed to be 
provided by electrical energy input to the feed gas mixture. The 
mixing chamber pressure was 4690 atm. A lower pressure, 2554 atm, 
was also tried for air and air + 10% NH 3 mixture. These values 
were chosen to keep the conditions used in this study similar to 
the conditions of the GASL study. In all computations reactions 
involving N0 2 were included and the rate data reported by 
Jachimowski (4) were used for these reactions. For all the other 
reactions GASL kinetics (1> were employed. 

The computation procedure was as follows: 

Find the initial mixture temperature to give the desired 
total enthalpy in the mixing chamber. Chemical and 
thermal equilibrium was assumed to exist in the mixing 
chamber. The equilibrium conditions and gas composition 
were found using the equilibrium combustion option of the 
GCKP . 

Pressure and temperature at the throat are calculated 
from the reservoir conditions found in the first step 
using 


P* =P 0 




) 


Gas composition at the throat was assumed to be the same 
as that in the mixing chamber. 

- The gas velocity was adjusted to give an initial Mach 
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number slightly above 1 (usually around 1.025). In some 
cases, due to the reactions taking place, just after the 
throat, temperature increases and Mach number decreases 
and the flow is choked not at the throat but at some 
distance from it in the expanding section of the nozzle. 
This necessitates starting Mach numbers even as high as 
1.16 at the throat. 

- Adjust the throat area to give the desired mass flow 
rate. 

Results: The results are shown in Figures 1 through 9. In each 

figure the first chart shows the variation of gas composition up to 
an axial distance of 1000 cm while the second chart covers the 
region up to an axial distance of 500 cm and is a blown up version 
of the first one. The results indicate that for the conditions 
specified, the flow in the nozzle is essentially frozen beyond 500 
cm and for this reason the flow conditions and gas compositions at 
500 cm are summarized in Table 2. 

The conclusions derived from Table 2 can be summarized as 
follows : 

- The addition of NH 3 or H 2 0 to air decreases the Mach number at 
the nozzle exit. This is mainly due to smaller y values of 
these mixtures. To achieve Mach 16 at the nozzle exit the 
expansion angle of the nozzle after about 300 cm from the 
throat should be increased gradually. An initial section with 
5 degree expansion is necessary to preserve the high rate of 
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Figure la 



chamber pressure 2554 atm. Air feed 





Figure lb 



amber pressure 






Figure 2a 



chamber pressure 2554 atm. 




Figure 2b 



Gas composition during reactive flow through 5 nozz 
chamber pressure 2554 atm. Feed is 10% NFL in air. 





Figure 3a 



chamber pressure 4690 atm. Air feed. 






ure 4690 atm. 


Figure 4a 



chamber pressure 4690 atm. 
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chamber pressure 4690 atm. 





Figure 5a 



chamber pressure 4690 atm. Feed is air + FLO equivalent to 




Figure 5b 



chamber pressure 4690 atm. Feed is air + H~0 equivalent to 





Figure 6a 



chamber pressure 4690 atm. Some N~ in air replaced to give 10% 
NH 0 in the feed. 




Figure 6b 



chamber pressure 4690 atm. Some 





Figure 7a 



chamber pressure 4690 atm. Air feed. 






Figure 8a 



chamber pressure 4690 atm. 





Figure 8b 



Gas composition during reactive flow through 10 nozzle 
chamber pressure 4690 atm. Feed is 10% NFL in air. 




Figure 9a 



chamber pressure 4690 atm. Feed is 10% NH, in oxygen-enriched air 




Figure 9b 



Gas composition during reactive flow through 5 nozzle. Mixing 
chamber pressure 4690 atm. Feed is 10% NH ? in oxygen-enriched air 





recombination reactions. 


Table 2. Summary of results for reactive nozzle flows for various 
gas compositions. Ma = Mach Number, y s = mole fraction 
of species i. 


Gas 

composi- 

tion 

Air 

Air + 10% NH 3 

Some N 2 in 
air replaced 
by NH 3 to 
give 10% NH 3 

air + H 2 0 
equivalent 
to H in 

10 % nh 3 

P 0 (atm) 

2254 

4690 

4690 

2254 

4690 

4690 

4690 

4690 

6 (deg .) 

5 

5 

10 

5 

5 

10 

5 

5 

Conditions at 500 cm from nozzle throat 

P (atm) 

.0076 

.0062 

.0009 

.0107 

.0087 

.0012 

.0086 

.0091 

T (K) 

636.5 

517.6 

320.8 

863.9 

703.4 

409.7 

714.2 

749.1 

Ma 

10.18 

11.37 

14.45 

8.89 

9.91 

12.95 

9.85 

9.68 

Y 

1.375 

1.386 

1.399 

1.337 

1.353 

1.384 

1.352 

1.346 

■■ 

.0174 

.0050 

.0070 

.0033 

.0013 

.0019 

.0014 

.0019 

YnO 

.0500 

.0535 

.0570 

.0333 

.0347 

.0381 

.0376 

.0434 


.1740 

.1797 

.1753 

.0953 

.0942 

.0932 

.1129 

.1578 

Mi 




.0086 

.0034 

.0059 

.0033 

.0034 

MB 




5.9E-4 

2.9E-4 

2.0E-4 

3.0E-4 

3.5E-4 

■Ml 




.0036 

.0018 

.0028 

.0016 

.0014 

MB 




.1363 

.1416 

.1357 

.1419 

.1385 

1 

.1000 

.0278 

.0399 

.0346 

.0138 

.0204 

.0124 

.0120 

y-o/Yo, 

.2874 

.2977 

.3252 

.3494 

.3684 

.4088 

.3330 

.2750 


— As expected, a larger nozzle expansion angle freezes the 
reactions earlier. As a consequence, the mole fractions of 0, 
NO, and H are larger while those of 0 2 and H 2 are smaller when 
the flow is expanded faster. The OH mole fraction is smaller 
for 10 degree expansion. The most significant changes are 
observed for the 0 mole fraction in pure air and H mole 
fraction in air containing 10% NH 3 . 

- When column 2 is compared to column 5 it is seen that an 
original gas mixture containing 10% NH 3 results in a test gas 
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with significantly lower O and NO mole fractions at the 
expense of having some H and OH in the test gas. The H mole 
fraction is over an order of magnitude smaller than that of NO 
and the OH mole fraction is over an order of magnitude smaller 
than that of H. Since the 0 2 mole fractions are different, to 
obtain a better comparison, 0/0 2 and N0/0 2 ratios are given in 
rows 15 and 16, respectively. These entries indicate that 
although the addition of ammonia to air is effective in 
reducing the 0/0 2 ratio, it increases the N0/0 2 ratio due to 
the presence of a larger number of N atoms. 

Since some air in the feed gas is replaced by NH 3 in the 10% 
ammonia containing feed, the 0 2 mole fraction in the test gas 
for this case is much lower and the number of N atoms present 
is much higher than with pure air feed. To. obtain a better 
comparison, the results for a feed gas in which some nitrogen 
is replaced by ammonia to give 10% NH 3 , are given in column 7. 
Even in this case the 0 2 mole fraction is lower than with pure 
air feed because some oxygen is used to form water. The 0/0 2 
ratio is lower than that obtained with pure air feed and the 
N0/0 2 ratio is only slightly higher. There was a significant 
improvement over the 10% NH 3 feed. The H and H 2 mole fractions 
are also lower with only insignificant increases in OH and 
H 2 0. Since the OH concentration continues to decrease up to 
the nozzle exit unlike the other species whose concentrations 
are frozen after about 500 cm after the nozzle throat, it is 
expected that the OH mole fraction in the test gas will be 
lower than the value reported in Table 2. 

To see if the beneficial effect of NH 3 is due to the H 2 0 
formed, the results obtained for a gas feed containing H 2 0 
equivalent to H in the 10% NH 3 feed are included in column 8 
of Table 2. These results show that the feed containing H 2 0 
gives the best 0/0 2 and N0/0 2 ratios. Even the H 2 and H 2 0 mole 
fractions are lower than those in column 6 with only slight 
increases in H and OH. 
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Table 3. Summary of results for reactive nozzle flows for various 
gas compositions. Mixing chamber pressure 4690 atm, 
nozzle expansion 5 degrees. Ma = Mach number, y ; = mole 
fraction of species i. 


Gas 

composi- 

tion 

Air 

Some N 2 in air replaced 
by NH 3 to give 10% NH 3 
+ .08 mol 0 2 per mol 

Air + H 2 0 
equivalent to H in 

10% nh 3 

P 0 (atm) 

4690 

4690 

4690 

0 (deg .) 

5 

5 

5 

Conditions at 500 cm from nozzle 

throat 

P (atm) 

.00618 

.00891 

.00911 

T (K) 

517.59 

732.05 

749.08 

Ma 

11.37 

9.759 

9.678 

Y 

1.3856 

1.3483 

1.3462 

Yo 

.004958 

.001811 

.001924 

y w 

.053517 

.045193 

.043374 

Yo 2 

.179664 

.171333 

.157830 

Yh 


.002975 

.003408 

Yoh 


.000324 

.000354 

Yh2 


.001204 

.001430 

Yh20 


.132001 

.138533 

Yo/Yo2 

.027596 

.010570 

.012190 

Y«/ Yoj 

.297873 

.263773 

.274815 


- To compare the results for feed gas compositions giving 
similar 0 2 mole fractions in the test gas, results were also 
obtained for a feed gas containing 10% NH 3 and enriched in 
oxygen to give an 0 2 mole fraction similar to that in the test 
gas obtained with pure air. These results are presented in 
column 2 of Table 3 and compared with some of the results from 
Table 2. These results indicate that better 0/0 2 and N0/0 2 
ratios can be obtained with oxygen enriched air containing 10% 
NH 3 than air containing water in an amount containing same 
number of H atoms as in the 10% NH 3 mixture. In addition, the 
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mole fractions of H, OH, H 2 , and H 2 0 are all lower with the 10% 
NH 3 containing enriched air. 


Conclusions 

The results shown above indicate that the addition of nitrogen 
and hydrogen-containing compounds to oxygen-enriched air can 
improve the 0/0 2 and N0/0 2 ratios and can have lower H, OH, H 2 , and 
H 2 0 mole fractions than with feed gases containing H 2 0. It appears 
that the liquid propellant HAN (N 2 H 4 0 4 ) may even be a better 
additive than NH 3 due to its lower N/H ratio and high oxygen 
content. Other propellants such as HMX, RMX, and HNTO will produce 
some C0 2 which will reduce the specific heat ratio (thus improving 
the required stagnation conditions) and increase the rates of 
recombination reactions in the nozzle. The presence of C0 2 in the 
test gas may partially counteract the effects of 0, H, and OH 
during combustion tests. 
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